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tion  with  Laser  Doppler  velocimetry  to  study  two  reacting  flows;  the  interaction 
of  a  Karman  vortex  street  with  a  flame,  and  the  interaction  of  grid-produced 
turbulence  with  a  flame.  Vortex  street  interaction  is  characterised  by  radical 
distortion  of  the  flame  front.  Dilatation  effects  downstream  of  the  flame 
front  dominate,  and  vortices  are  not  discernible  in  this  region.  Numerical 
modeling  gave  qualitatiive  agreement  with  the  experimental  results.  For  the 
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^compared  with  the  predictions  of  the  Bray-Moss-Libby  model.  It  was  found  that 
intermediate  states  had  to  be  taken  into  account.  Two-point  density  correlation 
measurements,  the  first  ever  reported#  showed  that  the  length  scales  of,  turbu¬ 
lence  within  the  reaction  cone  were  the  same  in  all  three  orthogonal  directions. 
Space-time  correlation  measurements  in  the  flow  direction  were  in  accord  with 
the  Taylor  hypothesis. 
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The  need  for  an  improved  understanding  of  the  effect  of  ~ 
turbulence  on  chemically  reacting  flows  has  long  been  recognized. 
However,  the  inherent  experimental  and  theoretical  difficulties 
associated  with  these  flows  have  continued  to  impede  progress 
towards  the  development  of  models  which  could  eventually  be 
used  in  practical  design  situations.  From  an  experimental 
point  of  view,  the  difficulty  of  obtaining  flow  measurements 
is  attested  to  by  noting  that  as  recently  as  1965  a  review  of 
experimental  data  (Williams,  1965)  showed  that  the  rms  velo¬ 
city  was  the  only  turbulent  parameter  reported.  New  measure¬ 
ment  techniques,  notably  laser  Doppler  velocimetry  (LDV)  have 
improved  the  situation,  but  the  experimental  data  available 
is  still  sparse  compared  with  data  for  other  turbulent  flows. 

The  theoretical  picture  for  turbulent  flows  with  combustion 
is  complicated  by  the  necessity  of  introducing  density  vari¬ 
ations  in  the  theoretical  specification  of  the  problem. 

Numerical  models  for  turbulent  flows  in  which  the  fluid  may 
be  considered  incompressible  now  successfully  predict  details 
of  fairly  complex  flows (Reynolds,  1976) .  However,  these 
models  are  not  easily  applied  to  combustion  problems  since 
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the  introduction  of  variable  density  results  in  terms  which 
have  no  counterpart  in  the  incompressible  models  (Libby  and 
Williams,  1976).  A  set  of  conservation  equations  whose  general 
form  is  similar  to  the  incompressible  time  averaged  equations 
results  if  density  weighted  time  averages  are  employed.  No 
data  exists  to  indicate  whether  analogous  terms  from  this  set 
of  conservation  equations  may  be  modelled  using  results  from 
incompressible  flows  as  a  guide.  In  order  to  investigate 
this  possibility  correlations  of  density  and  velocity  need  to 
be  measured. 

Under  previous  AFOSR  grants  feasibility  studies  indicated 
Rayleigh  scattering  to  be  a  reliable  and  accurate  method  of 
determining  gas  density  under  typical  combustion  conditions 
(Pitz  et  al.,  1976;  Cattolica  et  al.,  1977).  These  Rayleigh 
diagnostic  techniques  have  been  complemented  by  the  develop¬ 
ment  of  LDV  techniques  in  both  laminar  and  turbulent  flows. 

In  the  present  work  both  techniques  are  applied  to  the  study 
of  the  propagation  of  flame  fronts  in  turbulent  flows. 

Rayleigh  scattering  in  gases  results  from  the  inhomogeneous 
nature  of  the  medium  produced  by  fluctuations  in  the  dielectric 
constant.  The  Rayleigh  scattering  intensity  is  related  to  the 
gas  density  as  below: 


IR  =  C  IQ  N  Z  Xi  0R. 


where  IR  is  the  intensity  a 
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calibration  constant  of  the  optics,  IQ,  the  incident  laser 
light  intensity,  N,  the  total  molecular  number  density,  X^, 
the  mole  fraction  of  chemical  species,  and  aR^,  the  Rayleigh 
cross-section  for  the  ith  species.  Although  (1)  indicates  IR 
to  be  dependent  upon  the  degree  of  reaction  through  the 
species-dependent  terms ,  calculations  indicate  that  for  typical 
combustion  conditions,  the  Rayleigh  scattering,  to  an  accuracy 
of  better  than  3%,  is  independent  of  the  change  in  gas  compo¬ 
sition  and  directly  proportional  to  density.  This  is  primarily 
due  to  the  relatively  small  differences  in  crR^  and  the  dilute 
nature  of  products  and  fuel  components  in  air-fuel  mixtures. 

This  result  has  been  corroborated  by  comparing  theo¬ 
retical  results  with  measurements  in  the  laminar  boundary 
layer  of  a  heated  flat  plate  {Schefer  et  al.,  1980)  and  the 
successful  application  of  Rayleigh  scattering  to  laminar  flows 
with  heat  transfer  and  combustion  has  led  naturally  to  the 
initiation  of  studies  of  disturbed  and  turbulent  combustion 
flows. 

During  the  period  of  the  subject  AFOSR  contract,  combined 
LDV  and  Rayleigh  diagnostic  techniques  were  used  to  study  the 
following  problems. 

(1)  The  interaction  of  a  K£rm5n  vortex  street  with  a 
flame  front; 

(2)  The  interaction  of  grid-produced  turbulence  with  a 


flame  front. 


.  The  interaction  of  a  K^rm^n  vortex  street  with  flame  front 


The  interest  in  this  problem  derives  from  the  observation 
that  large-scale  'coherent'  structures  play  an  important  role 
in  determining  the  characteristics  of  flame  propagation  in 
turbulent  flows.  The  KSrm^n  vortex  street  provides  an  ideal¬ 
ization  of  this  interaction.  The  results  of  this  study  are 
reported  in  Namer  (1980)  and  Namer  et  al.  (1982).  They  may 
be  summarized  as  follows: 

a)  The  mean  statistical  properties  of  the  flame  show  the 
development  of  the  turbulent  flame  brush  by  the  vortex  street 
and  the  reduction  in  turbulent  intensity  behind  the  flame 
along  with  the  disappearance  of  a  discernable  wake  from  the 
turbulence-generating  rod.  The  probability  density  functions 
of  the  density  indicate  the  general  validity  of  a  flame  sheet 
model  based  on  the  presence  of  burned  and  unburned  gas,  but 
with  significant  (~  30%)  contribution  of  intermediate  states. 
Detailed  comparison  with  some  of  the  theoretical  predictions 
of  Bray  and  Libby  (1976)  and  Bray  and  Moss  (1977)  show  reason¬ 
able  agreement  but  also  indicate  the  necessity  of  taking 
intermediate  states  into  account  in  order  to  obtain  quanti¬ 
tative  agreement. 

b)  By  using  a  phase-locked  technique  for  recording  the  vortex 
structure  it  was  possible  to  obtain  complete  mapping  of  the 
flow  field  for  given  positions  of  the  KSrmcin  vortices  shed  by 
the  rod.  Contour  plots  of  constant  velocity  contours 
identified  the  vortex  structure.  When  a  flame  is  present  the 


vortex  structure  is  distorted  and  displaced  upstream  of  the 
flame  front,  and  is  not  discernable  downstream  of  the  flame. 
The  flame  front  is  radically  distorted  by  the  vortices  as  is 
most  clearly  evident  from  the  density  measurements.  As 
expected  from  the  scale  of  the  'turbulence1  and  the  Darakohler 
number,  it  appears  that  there  is  only  a  single  flame  front 
which  is  distorted,  and  that  pockets  of  burned  and  unburned 
gas  are  not  being  generated. 

The  experimental  results  were  compared  with  numerical 
calculations  of  the  flow  and  flame  front  due  to  Karasalo  et 
al.  (1980).  In  these  calculations  single  vortices  are 
inserted  into  the  flow  to  simulate  the  K^rmSn  vortex  street, 
and  the  flame  is  modelled  as  a  line  source  of  specific  volume. 
There  was  qualitative  agreement  between  the  experiment  and 
the  model  calculations,  but  the  calculations  indicated  con¬ 
siderably  less  perturbation  of  the  flame  front  than  found 
experimentally . 

2.  The  interaction  of  grid-produced  turbulence  with 
a  flame  front 

An  extensive  set  of  measurements  on  turbulent  flame 
structure  in  the  presence  of  homogeneous,  isotropic  grid¬ 
generated  turbulence  have  been  carried  out.  The  results  are 
reported  in  Bill  et  al,  (1981),  Bill  et  al.  (1982)  and 
Namazian  et  al.  (1982) 

In  the  first  set  of  measurements,  LDV  data  were  combined 


with  Rayleigh  data  to  obtain  density  and  streamwise  velocity 
statistics  for  a  pre-mixed  V-shaped  ethylene/air  flame.  Tur¬ 
bulent  flame  propagation  was  studied  for  three  upstream  flow 
velocities,  U  =  2.45,  5.0  and  6.84  m/s,  with  equivalence 
ratios  ranging  from  0.55  to  0.75.  The  angle  of  the  flame  with 
the  flow  axis  was  varied  from  12°  to  24°.  Profiles  of  the 
mean  values  and  fluctuation  intensities  of  velocity  and  density 
were  obtained  through  the  flame.  For  all  conditions,  the  velo¬ 
city  fluctuation  intensity  was  observed  to  decrease  behind  the 
flame,  indicating  the  dominance  of  dilatation  effects.  Addi¬ 
tionally,  the  profiles  of  the  mean  streamwise  velocity  component 
indicated  that  the  flow  was  significantly  deflected  at  the  flame 
front,  particularly  for  the  case  of  the  most  oblique  (12°)  flame. 
In  this  flame  condition  the  wake  of  the  flameholder  was  found 
to  dominate  the  flowfield  behind  the  flame  and  vortices  were 
observed  in  the  burned  gas.  Thus  the  results  for  this  flame 
configuration  were  not  typical  of  the  interaction  of  a  flame 
with  grid-produced  turbulence.  Probability  density  functions  of 
the  Rayleigh-scattered  radiation  and  measured  turbulence  inten¬ 
sities  indicated  that  states  intermediate  between  the  unburned 
and  burned  gases  are  likely  to  be  significant  in  determining 
density  statistics. 

In  the  second  set  of  measurements  a  more  detailed  comparison 
is  made  between  the  Rayleigh  scattering  results  for  the  density 
statistics  in  a  premixed  V-shaped  flame  and  the  predictions  of 
the  Bray-Moss-Libby  model.  Probability  density  distributions 
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of  density  fluctuations  are  obtained  at  a  number  of  locations 
through  the  flame.  A  comparison  of  these  results  with  the 
predictions  of  the  B-M-L  model,  in  which  intermediate  states 
are  neglected,  shows  that  this  neglect  is  not  justified  for 
the  conditions  of  the  experiment.  However',  the  measured 
skewness  of  the  density  fluctuations  is  accurately  predicted 
by  the  model.  The  data  suggest  that  it  should  be  possible 
with  relatively  modest  modifications  in  the  model  to  obtain 
satisfactory  agreement  between  theory  and  experiment. 

In  the  most  recent  set  of  measurements  on  V-shaped  flames 
propagating  into  grid-produced  turbulent  flow,  two- point 
correlations  of  density  have  been  obtained.  These  results 
represent  an  entirely  new  advance  in  the  diagnostics  of  turbu¬ 
lent  flows.  Since  they  have  not  yet  been  published,  we  report 
them  in  detail. 

A  unique  beam  splitter  plus  split-lens  optical  system  was 
employed  to  produce  two  laser  beams  whose  waist  locations 
could  be  varied  in  a  continuous  manner  from  coincidence  to  any 
arbitrary  distance  apart  (limited  only  by  the  diameter  of  the 
optics) .  Separate  photomultipliers  focused  on  the  waists  of 
each  of  the  beam  simultaneously  recorded  the  Rayleigh 
scattering  at  the  waist  locations.  These  signals  were  sampled 
at  a  10  kHz  rate  with  a  PDP-11/10  computer  and  stored  on 
magnetic  tape  for  post-processing.  The  flow  system  is  shown 
in  Figure  1  and  the  optical  system  is  shown  in  Figure  2. 

The  two  point  correlation  function  R(«S)  is  defined  by 


R(6) 


<p^ (t) (t) > 
<pJ(t)2>1/2<P2(t) 


57T7? 


where  6  is  the  distance  of  separation  of  the  measurement 
positions  denoted  by  subscripts  1  and  2,  and  p'  is  the  density 
fluctuation  and  <>  denotes  the  time  average.  Measurements  of 
R( 6)  were  made  through  the  flame  with  6  taken  in  each  of  the 
three  coordinate  directions,  x,  y,  z  (Figure  1).  The  results 
are  shown  in  Figure  3.  It  is  seen  that  the  values  of  the 
spatial  correlations  are  very  similar  in  the  three  directions, 
indicating  that  the  spatial  structure  of  the  turbulent  flame, 
under  the  conditions  of  the  experiment,  is  such  that  the 
associated  integral  scales  are  the  same  in  these  directions. 

The  integral  scale,  represented  by  the  area  under  the  R(<5) 
curve,  is  found  to  be  approximately  2  mm,  about  half  the  thick¬ 
ness  of  the  turbulent  flame  zone.  Interestingly,  the  trans¬ 
verse  integral  scale  of  the  grid  turbulence  was  also  about 
2  mm. 

The  space-time  correlation,  defined  by 


R(5,t)  = 


<PX (t) 


<P^(t)  p^(t  +  t)> 

27T72 


2^757  . 


C2(t) 


where  x  is  the  time  delay,  was  also  evaluated,  and  is  shown 
in  Figure  4,  as  a  function  of  x  for  several  values  of  6  taken 
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along  the  x  or  flow  direction.  The  maxima  in  each  curve  is 
proportional  to  6,  with  the  proportionality  factor  about  7  m/s, 
which  was  the  flow  velocity.  Since  previous  measurements  have 
indicated  that  for  the  conditions  of  the  experiment  the  axial 
velocity  does  not  change  significantly  through  the  flame, 
this  indicates  that  the  turbulent  density  structures  are  simply 
convected  downstream  by  the  flow,  as  is  suggested  by  the 
Taylor  hypothesis. 

In  Figure  5  the  maximum  value  of  R(6,x)  is  plotted  as  a 
function  of  6.  The  decrease  of  R(5,x)  with  6  implies  that  the 
turbulent  density  structures  decay  as  they  are  convected  by 
the  flow.  Figure  6  shows  R(6,x)  for  6  taken  in  the  transverse 
direction  within  the  flame.  In  this  case  the  maximum  is 
always  found  at  x  =  0. 

All  of  these  results  should  be  of  use  in  improved  modelling 


of  turbulent  flows  with  combustion 


10 


1 


REFERENCES 

Bill,  R.  G.,  Jr.,  Namer,  I.  and  Talbot,  L.  1981,  Combust. 
Flame  4_ 3,  229. 

Bill,  R.  G.,  Jr.,  Namer,  I.,  Talbot,  L.  and  Robben,  F.  1982, 
Combust.  Flame  4^,  277. 

Bray,  K.N.C.  and  Libby,  P.  A.  1976,  Phys.  Fluids  19,  1687. 
Bray,  K.N.C.  and  Moss,  J.  B.  1977  Acta  Astronautica  4,  291. 


Cattolica,  R. ,  Robben,  F  and  Talbot,  L.  1977  in  "Experimental 
Diagnostics  in  Gas  Phase  Combustion  Systems,"  Prog,  in 
Astr.  and  Aero.  5_3,  AIAA ,  New  York. 

Karasalo,  I.,  Chorin,  A.  J.  ,  Namer,  I.,  Robben,  F.  and  Talbot,  L. 
1980,  Lawrence  Berkeley  Laboratory  Report  L3L-10S79. 

Libby,  P.  A.  and  Williams,  F.  A,  1976  in  Annual  Reviews  of 
Fluid  Mechanics  8. 

Namazian,  M. ,  Talbot,  L.  and  Robben,  F.  1982.  "Two  Point 

Rayleigh  Scattering  Measurements  in  a  V-shaped  Turbulent 
Flame."  Submitted  to  the  19th  International  Symposium 
on  Combustion,  Haifa,  Israel. 

Namer,  I.  1980,  "An  Experimental  Investigation  of  the  Inter¬ 
action  between  a  Karmdn  Vortex  Street  and  a  Premixed 
Flame,"  Ph.D.  Thesis,  University  of  California,  Berkeley. 

Namer,  I.,  Bill,  R.  G.,  Jr.,  Talbot,  L.  and  Robben,  F.  1982, 

AIAA  J.  (to  appear) . 

Pitz,  R.  W.  ,  Cattolica,  R.  ,  Robben,  F,  and  Talbot,  I..  1976, 

Combust.  Flame  2J_,  313. 

Reynolds,  W.  C.  1976,  "Computation  of  Turbulent  Flows"  in 

Annual  Reviews  of  Fluid  Mechanics  j3,  Annual  Reviews,  Inc., 
Palo  Alto. 

Schefer,  R.  W. ,  Robben,  F.  and  Cheng,  R.  K.  1980,  Combust. 

Flame  38,  51. 


Williams,  F.  A.  1965,  "Combustion  Theory,"  Addison-Wesley 
Reading,  Mass. 


11 


Fig,  1.  Schematic  of  the  experimental  apparatus  for 
two  point  Rayleigh  scattering  measurements. 


Fig.  2.  The  optical  set-up  for  two  point  Rayleigh  scattering  technique.  Single  beam 

SYSTEM  PROVIDES  SEPARATION  BETWEEN  POINTS  IN  HORIZONTAL  DIRECTION  AND  DUAL  BEAM  SYSTEM 
IN  VERTICAL  DIRECTION. 
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DISTANCE  BETWEEN  SAMPLED  POINTS  (mm),  & 

Fig.  3.  The  density  spatial  correlation  in  a  turbulent  flame  in  three  orthogonal  X/ 

Y,  AND  Z  DIRECTIONS.  MEASUREMENT  POINT  35  MM  DOWNSTREAM  FROM  THE  FLAME  HOLDER. 


THE  DISTANCE  BETWEEN  SAMPLED  POINTS. 
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Time^space  correlation  versus  the  delay  time  for  separation 


